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Project  Abstract 

The  goal  of  this  proposal  is  to  model  the  stability  of  quantum  states  of  matter  derived 
from  topological  insulators  against  two  types  of  corrections:  strong  inter-particle  inter¬ 
actions  and  heating.  I  will  examine  interacting  atoms  in  square  optical  lattices  with  spin 
orbit  coupling,  and  more  generally,  gauge  fields,  as  a  route  to  building  Hubbard  models 
hosting  fractional  topological  insulators.  I  will  analyze  these  models  by  combining  nu¬ 
merical  exact  diagonalization  on  small  lattice  clusters  with  an  analytic  variational  theory. 
I  also  propose  to  study  a  new  model  of  finite  temperature  topological  superconductors 
of  dipoles  placed  in  an  optical  lattice.  I  will  construct  and  analyze  a  model  using  a  com¬ 
bination  of  mean  field  theory  and  quantum  Monte  Carlo.  The  proposed  work  will  foster 
new  directions  in  experiments  with  optical  lattices  containing  cold  atomic  gases.  These 
new  states  of  matter  should  exhibit  new  particles  as  excitations.  Analyses  of  the  stability 
of  these  new  topological  phases  will  thus  play  a  crucial  role  in  advancing  fundamentally 
new  directions  in  physics. 
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(1)  List  of  Appendices,  Illustrations  and  Tables 


•  Figure  1:  Schematic  of  dipolar  fermions  (spheres)  in  a  2D  optical  lattice.  Dipolar 
moments  (arrows  on  each  sphere)  align  along  an  applied  field,  at  an  angle  9  with 
the  x-axis.  The  optical  lattice  depth  is  large  along  the  y  direction. 

•  Figure  2:  The  thermal  expectation  value  of  string  operators  from  quantum  Monte 
Carlo  as  a  function  of  an  applied  global  field  for  several  system  sizes.  The  top  (bot¬ 
tom)  panel  shows  data  for  a  characteristic  low  (high)  temperature 

•  Figure  3:  The  density  of  a  "vortex  core"  ansatz  state  of  helical  fermions  in  the  pres¬ 
ence  of  Rashba  spin  orbit  coupling  plotted  as  a  function  of  x  and  y  in  units  of  the 
spin-orbit  length  /so. 

•  Figure  4:  Left:  Plot  of  the  many-body  energy  of  a  diagonal  nearest  neighbor  inter¬ 
action  as  a  function  of  wavevector.  The  ground  state  is  a  Wigner  crystal  and  the 
excitations  are  classical.  Right:  The  same  but  for  a  flat  spin-orbit  band.  The  ground 
state  is  still  a  Wigner  crystal  but  the  excitations  are  show  an  emergent  dispersion 
due  to  quantum  effects. 

•  Figure  5:  Spin  structure  factors  for  different  quantum  phases.  The  upper  panels 
show  the  results  for  fermions,  while  the  lower  panels  show  the  results  for  bosons. 


(2)  Statement  of  The  Problem  Studied 

The  goal  of  the  STIR  proposal  was  to  study  the  stability  of  quantum  states  of  matter 
derived  from  topological  insulators  against  two  types  of  corrections:  strong  inter-particle 
interactions  and  heating.  These  parameters  can  be  tuned  and  controlled  in  optical  lattices 
while  leaving  out  other  influences  by  construction. 

We  studied  strong  spin  orbit  coupled  problems  to  set  the  stage  for  further  work  on 
topological  states  of  matter.  Topological  states  of  matter  are  derived  from  strong  spin  or¬ 
bit  interaction  in  the  context  of  topological  insulators  and  topological  superconductors. 
In  section  3A  I  discuss  work  that  studies  the  thermal  stability  of  topological  supercon¬ 
ductors,  in  section  3B  I  discuss  work  that  builds  a  formalism  to  study  the  competition 
between  strong  interaction  and  strong  spin  orbit  coupling.  In  section  3C  I  discuss  results 
on  the  search  for  topological  states  in  the  presence  of  strong  spin  orbit  coupling  and  strong 
interactions  between  fermions  in  an  optical  lattice.  And  finally,  in  section  3D,  I  discuss 
results  that  studies  strong  interaction  and  strong  spin  orbit  coupling  between  bosons  in 
an  optical  lattice. 

The  work  started  here  can  be  used  to  study  topological  states  of  matter  in  more  re¬ 
alistic  settings.  Comparison  with  and  predictions  for  experiments  are  now  a  possibility. 
Furthermore  the  methods  and  techniques  constructed  here  can  be  generalized  to  apply  to 
other  systems:  quantum  wire  arrays  containing  topological  superconductors,  quantum 
Hall  effects,  fractional  Chern  insulator  lattice  models,  and  fractional  topological  insula¬ 
tors. 


(3)  Summary 

3A:  Enhancing  the  thermal  stability  of  entanglement  between  Major  ana  fermions  with  dipoles  in 
optical  lattices 


The  wide  variety  of  optical  lattice  ge¬ 
ometries  offer  unprecedented  tunability  in 
manipulating  quantum  degenerate  gases 
into  complex  quantum  states.  Recent  de¬ 
velopments  in  the  cooling  of  molecules 
and  magnetic  atoms  imply  that  anisotropy 
in  dipolar  interactions  will  soon  provide 
further  opportunity  to  explore  some  of 
the  most  elusive  yet  compelling  quan¬ 
tum  states,  entangled  Major  ana  fermions 
(MFs).  Fig.  1  shows  a  figure  of  a  lattice 
model  we  studied. 

Certain  models  demonstrate  particle¬ 
like  excitations  that  behave  as  MFs  thanks 
to  peculiar  non-local  symmetries.  MFs 
are  their  own  antiparticles.  They  en¬ 
tangle  with  each  other  over  large  dis¬ 
tances  through  string  operator  correla¬ 
tions,  signaling  underlying  topological  or¬ 
der  with  fascinating  properties  that  have 
motivated  proposals  for  topologically  pro¬ 
tected  qubits.  The  crossing  of  string  oper¬ 
ators  is  responsible  for  unusual  anyonic  braid  statistics.  And  string  operators  connecting 
these  excitations  also  underly  theories  of  quantum  state  teleportation. 

The  zero-temperature  properties  of  models  hosting  topological  order  set  the  stage  for 
work  connected  to  experiments.  Kitaev's  two-dimensional  (2D)  Toric  Code  Flamiltonian 
motivated  early  proposals  in  optical  lattices  containing  ultracold  atoms,  polar  molecules, 
and  atoms  in  Rydberg  states.  But  the  ID  Kitaev  chain  model  is  one  of  the  simplest  mod¬ 
els  supporting  free  MF  excitations.  Anticipation  of  non-local  MF  properties  in  ID  led 
to  experimental  proposals  in  both  optical  lattices  and  solids,  and  recent  experiments  in 
quantum  wires.  But  prospects  for  observing  interesting  effects  due  to  non-local  correla¬ 
tion  of  MF  pairs  over  long  times  and  distances  hinge  on  the  stability  of  string  operators. 

String  operators  in  important  lattice  models  are  unstable  at  non-zero  temperatures, 
e.g.,  string  operators  in  the  2D  Toric  Code  model  vanish  at  long  times  and  distances 
because  of  thermal  excitations.  Recent  work  also  argues  that  MFs  in  lattice  models  of 
topological  p-wave  superconductors  are  sensitive  to  thermal  fluctuations. 

We  addressed  the  issue  of  thermal  stability  by  showing  that  anisotropic  interactions 
between  dipolar  fermions  in  optical  lattices  can  be  used  to  significantly  enhance  ther¬ 
mal  stability  [1].  We  constructed  a  model  of  oriented  dipolar  fermions  in  a  square  opti¬ 
cal  lattice.  Our  model  reduces  to  an  interacting  Kitaev  chain  model  of  fermions  hoping 
with  an  energy  t  along  chains.  We  found  that  domains  established  by  strong  interactions 


Figure  1:  Schematic  of  dipolar  fermions 
(spheres)  in  a  2D  optical  lattice.  Dipolar  mo¬ 
ments  (arrows  on  each  sphere)  align  along 
an  applied  field,  at  an  angle  6  with  the  x-axis. 
The  optical  lattice  depth  is  large  along  the  y 
direction. 


exhibit  enhanced  correlation  between  Major  ana  fermions  over  large  distances  and  long 
times  even  at  finite  temperatures,  suitable  for  stable  redundancy  encoding  of  quantum 
information.  Our  approach  can  be  generalized  to  a  variety  of  configurations  and  other 
systems,  such  as  quantum  wire  arrays. 

Fig.  2  shows  one  of  our  key  results.  The 
expectation  value  of  the  string  operators, 

Pj,  act  as  order  parameters.  Unique  val¬ 
ues,  (Pj)  =  ±1,  can  be  used  to  define 
each  sector  and  therefore  indicate  stabil¬ 
ity  in  the  string  operators  of  MFs.  But 
(P)  =  0  indicates  that  thermal  excita¬ 

tions  destroy  any  distinction  between  sec¬ 
tors.  We  compute  (Pj)  explicitly  on  a 
model  describing  strongly  interacting  Ki- 
taev  wires  to  show  indications  of  sponta¬ 
neous  breaking  of  these  discrete  symme¬ 
tries  for  even  at  non-zero  temperatures. 

To  detect  such  a  symmetry  breaking  we 
perturb  the  above  spinless  fermion  model 
with  a  weak  global  field  of  strength  h.  The 
global  field  imposes  a  splitting  between 
the  otherwise  degenerate  states.  Our  re¬ 
sults  clearly  show  the  thermal  robustness 
of  Majorana  fermions  in  our  model. 

We  considered  an  effective  model  of 
oriented  dipolar  fermions  in  a  2D  lat¬ 
tice  that  allows  hopping  along  directions 

where  the  dipoles  attract  but  suppresses  hopping  along  directions  where  dipoles  repel. 
We  also  performed  unbiased  QMC  to  directly  check  that  string  operators  defining  non¬ 
local  MF  states  remain  robust  to  thermal  fluctuations.  This  class  of  models  also  applies  to 
Coulomb-coupling  in  MF  models  of  quantum  wire  arrays  or  quasi-ID  tubes  containing 
topological  superconductors. 


Figure  2:  The  thermal  expectation  value  of 
string  operators  from  quantum  Monte  Carlo 
as  a  function  of  an  applied  global  field  for 
several  system  sizes,  L.  The  top  (bottom) 
panel  shows  data  for  a  characteristic  low 
(high)  temperature. 


3B:  Vortex  Attachment  via  Matrix  Products:  Application  to  Flat  Spin-Orbit  Bands 

Strongly  correlated  systems  can  often  be  understood  using  wavefunctions  constructed 
from  Jastrow  factors  that  lower  interaction  energies  by  separating  particles.  A  useful 
interpretation  of  Jastrow  factors  based  on  vortices  (wavefunction  zeroes)  views  a  Jastrow 
factor  as  tool  to  introduce  vortices  into  an  otherwise  weakly  correlated  system.  Careful 
tuning  of  the  location  of  each  vortex  optimizes  energetics  of  individual  problems.  This 
procedure  has  proven  to  be  quantitatively  useful  in  combination  with  stochastic  methods 
in  quantum  chemistry  and  other  fields. 

Vortex  insertion  with  Jastrow  factors  has  been  the  tool  of  choice  in  highly  successful 
theories  of  the  quantum  Hall  regime.  In  the  quantum  Hall  regime  the  kinetic  energy 
is  quenched  leaving  the  interaction  to  operate  in  a  flat  band.  The  Hamiltonian  has  no 
small  parameters  and  thus  defines  a  challenging  many-body  problem.  In  this  regime 
the  success  of  quantum  Hall  ansatz  wavefunctions  shows  that  vortex  attachment  to  each 
particle  helps  screen  the  otherwise  strong  inter-particle  repulsion  by  tuning  the  size  of 
correlation  holes. 

Flat  band  models  have  recently  been 
shown  to  describe  a  wealth  of  other  quan¬ 
tum  many-body  systems.  These  include 
graphene-based  nanostructures,  fast  rotat¬ 
ing  atomic  gases,  kagome  lattices,  frac¬ 
tional  Chern  insulators,  and  strongly  spin- 
orbit  coupled  systems,  to  name  a  few.  The 
success  of  the  quantum  Hall  wavefunc¬ 
tions  in  capturing  the  essential  physics  of 
the  Coulomb  interaction  within  a  flat  Lan¬ 
dau  level  suggests  that  vortex  attachment 
to  particles  might  be  a  quantitatively  use¬ 
ful  ansatz  in  solving  these  other  flat-band 
problems  as  well.  Vortex  attachment  may 
also  be  used  to  study  fractional  topological 
insulators. 

We  constructed  a  formalism  that  uses 
a  basis-independent  representation  of  Jas¬ 
trow  factors  to  attach  vortices  in  the  matrix  product  setting  [2],  We  then  tested  this  formal¬ 
ism  on  a  few  examples.  We  showed  that  the  formalism  can  be  used  to  directly  write  the 
single  component  Laughlin  states  in  terms  of  matrix  products.  We  also  used  the  formal¬ 
ism  to  solve  a  flat-band  problem  that  does  not  have  a  polynomial  basis  representation: 
two-dimensional  helical  fermions  in  the  presence  of  strong  Rashba  spin-orbit  coupling. 
This  problem  is  of  direct  relevance  to  ongoing  experiments  with  ultracold  atomic  gases. 
We  found  that  interactions  in  a  flat  spin-orbit  band  favor  the  formation  of  a  central  vortex 
in  analogy  to  what  has  been  found  in  studies  of  rotating  quantum  gases  and  quantum 
dots.  The  formalism  allows  straightforward  validation  of  Jastrow-correlated  ansatz  states 
in  small  system  sizes  and  sets  the  stage  for  optimization  of  these  wavefunctions  with  effi¬ 
cient  matrix  product  algorithms.  We  are  currently  preparing  a  manuscript  for  submission 
[2], 


Figure  3:  The  density  of  a  "vortex  core" 
ansatz  state  of  helical  fermions  in  the  pres¬ 
ence  of  Rashba  spin  orbit  coupling  plotted  as 
a  function  of  x  and  y  in  units  of  the  spin-orbit 
length  Iso ■ 


The  primary  result  was  a  comparison  of  ansatz  wavefunctions  and  results  from  exact 
diagonalization.  The  overlaps  between  the  ansatz  and  exact  states  were  better  than  98%. 
Fig.  3  shows  the  density  in  the  plane  for  the  successful  state.  The  vortex  in  the  system 
center  helped  minimize  the  interaction  energy. 


3C:  Flat-band  engineering  of  interactions  in  spin-orbit  coupled  ID  optical  lattices 


Emergent  excitations  can  have  properties  that  are  not  restricted  by  common  symme¬ 
tries.  One  of  the  best  known  examples  is  in  the  fractional  quantum  Hall  regime  where 
many-body  effects  in  a  quantum  flat-band  lead  to  excitations  with  topological  properties. 
These  topological  excitations  can  have,  for  example,  fractional  charge. 

Recent  cold  atom  experiments  have 
made  an  important  breakthrough  in  real¬ 
izing  spin-orbit  couplings,  which  opened 
a  door  to  realizing  flat  quantum  bands. 
We  find  that  strong  interactions  in  this  flat 
band  leads  to  intriguing  excitations  with 
fractional  charge. 

We  are  currently  preparing  our  manusc¬ 
ript  to  report  results  that  explores  these 
issues  in  optical  lattices  [3].  We  con¬ 
structed  a  tight-binding  Hamiltonian  for 
a  spin-orbit  coupled  Fermionic  optical  lat¬ 
tice  system,  and  perform  flat-band  pro¬ 
jection  of  on-site  interactions  to  gener¬ 
ate  longer-range  interactions  between  flat- 
band  particles.  We  solved  the  resultant 
projected  interacting  Hamiltonian  with  ex- 
Figure  4:  Feft:  Plot  of  the  many-body  energy  ac*  diagonalization  at  some  parameter 
of  a  diagonal  nearest  neighbor  interaction  as  re§ime'  and  obtaincd  sPinor  Wi§ner  crys‘ 


a  function  of  wavevector.  The  ground  state 
is  a  Wigner  crystal  and  the  excitations  are 
classical.  Right:  The  same  but  for  a  flat  spin- 
orbit  band.  The  ground  state  is  still  a  Wigner 
crystal  but  the  excitations  are  show  an  emer¬ 
gent  dispersion  due  to  quantum  effects. 


tal  ground  state  with  an  emergent  ki¬ 
netic  bands  above  the  ground  state.  We 
deduce  an  effective  extended  Hubbard 
model  with  modulated  next  nearest  neigh¬ 
bor  hoppings  to  describe  the  emergent  ki¬ 
netic  bands.  Emergent  dispersive  bands 
are  shown  in  Fig.  4.  The  effective  model 


then  leads  to  charge  fractionalization  of  the 
elementary  Bosonic  excitations,  revealed  in  the  conventional  bosonization  approach. 


3D:  Dzyaloshinskii-Moriya  Interaction  and  Spiral  Order  in  Spin-orbit  Coupled  Optical  Lattices 
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In  this  work  [4]  we  showed  that  the  recent  experimental  realization  of  spin-orbit  cou¬ 
pling  in  ultracold  atomic  gases  can  be  used  to  study  different  types  of  spiral  order  and 
resulting  multiferroic  effects.  Spin-orbit  coupling  in  optical  lattices  can  give  rise  to  the 
Dzyaloshinskii-Moriya  spin  interaction  which  is  essential  for  spin  spiral  order.  We  de¬ 
rived  an  effective  spin  model  in  the  deep  Mott  insulator  region  at  half  filling,  and  demon¬ 
strate  that  the  Dzyaloshinskii-Moriya  interaction  in  optical  lattices  can  be  made  extremely 
strong  with  realistic  experimental  parameters.  The  rich  phase  diagrams  of  the  effective 
spin  model  for  fermion  and  bosons  are  obtained  via  classical  Monte  Carlo  simulations. 

The  interplay  between  ferro¬ 
electric  and  ferromagnetic  order  5=0.13 

in  complex  multiferroic  materi¬ 
als  presents  a  set  of  compelling  ^ 
fundamental  condensed  matter 
physics  problems  with  potential 
multifunctional  device  applica¬ 
tions.  Ferroelectric  and  ferromag-  ^ 
netic  order  compete  and  normally 
cannot  exist  simultaneously  in 
conventional  materials.  Coexis¬ 
tence  of  these  orders  often  relies 
on  strong  correlation.  In  some  R  5.  S pin  stmcture  factors  for  different  quan- 

strongly  correlated  systems,  such  tum  phases  The  upper  panels  show  the  results  for 

as  the  perovskite  transition  metal  c  '  T  ,  ,u  c 

.  ,  r  ,  ,  termions,  while  the  lower  panels  show  the  results  tor 

oxides,  these  two  phenomena  can  . 

,,  1  .  bosons, 

occur  simultaneously.  Construc¬ 
tion  and  design  of  high- Tc  mag¬ 
netic  ferroelectrics  is  still  an  open  and  active  area  of  research.  These  and  other  complex 
solids  incorporate  different  types  of  interactions,  including  electron-electron  interactions, 
electron-phonon  interactions,  spin-orbit  couplings,  lattice  defects,  and  disorder,  making 
the  determination  of  multiferroic  mechanisms  a  remarkable  challenge  for  most  materi¬ 
als.  Due  to  the  inherent  complexity  of  the  different  interactions  involved,  the  underlying 
mechanism  for  the  emergence  of  multiferroic  behavior  is  not  fully  clear.  An  unbiased 
and  direct  method  to  explore  multiferroic  behavior  (e.g.,  dynamics,  critical  exponents, 
the  impact  of  controlled  disorder,  etc.)  in  an  ideal  setting  is  thus  highly  appealing. 

We  showed  that  the  power  of  optical  lattice  systems  to  emulate  magnetism  can  be 
combined  with  recent  experimental  developments  realizing  spin-orbit  coupling  to  emu¬ 
late  multiferroic  behavior.  Recently,  spin-orbit  coupled  optical  lattices  have  been  shown 
to  possess  flat  energy  band,  and  important  experimental  progress  has  been  made  for  ob¬ 
serving  the  single  particle  band  dispersion  (including  the  flat  band)  in  spin-orbit  coupled 
optical  lattices.  The  main  findings  of  this  work  are  the  following:  (I)  We  derive  the  effec¬ 
tive  Hamiltonian  for  spin-1 /2  fermions  and  bosons  in  optical  lattices  in  the  large  interac¬ 
tion  limit.  We  show  that  spin-orbit  coupling  leads  to  an  effective  in-plane  Dzyaloshinskii- 
Moriya  term,  an  essential  ingredient  in  models  of  spiral  order  and  multiferroic  effects  in 
general.  We  estimate  that  the  DM  term  is  of  the  same  order  as  the  Heisenberg  coupling 


constant,  /.  (II)  We  study  the  phase  diagram  of  the  effective  spin  model  using  classical 
Monte  Carlo.  We  find  that  competing  types  of  spiral  order  depend  strongly  on  the  spin- 
orbit  coupling  strength  and  effective  Zeeman  field.  Fig.  5  shows  examples  of  the  spin 
structure  factor  used  to  identify  phase  transitions.  (Ill)  We  find  that  the  critical  temper¬ 
ature  for  the  spiral  order  can  be  made  rather  high.  Thus,  if  magnetic  quantum  phase 
transitions  can  be  emulated  in  optical  lattices,  then  spiral  order  and  multiferroics  can  also 
be  realized  in  the  same  setup  with  the  inclusion  of  spin-orbit  coupling. 
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